Neutrino Bremsstrahlung in Neutron Matter from Effective Nuclear Interactions 
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We revisit the emissivity from neutrino pair bremsstrahlung in neutron-neutron scattering, 
nn nn uv, which was calculated from the one-pion exchange potential including correlation effects 
by Friman and Maxwell. Starting from the free-space low-momentum nucleon-nucleon interaction 
Viowfe, we include tensor, spin-orbit and second-order medium-induced non-central contributions 
to the scattering amplitude in neutron matter. We find that the screening of the nucleon-nucleon 
interaction reduces the emissivity from neutrino bremsstrahlung for densities below nuclear mat- 
ter density. We discuss the implications of medium modifications for the cooling of neutron stars 
via neutrino emission, taking into account recent results for the polarization effects on neutron 
superfluidity. 
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" The observation of neutron star properties and their evolution provides a challenging astrophysical setting for 
the study of dense nuclear matter [1,2]. On the theoretical side, the main objectives are a more comprehensive 
P\J understanding of the equation of state and of microscopic nuclear properties, such as pairing and transport phenomena, 
J> \ where progress is tied to improving many-body calculations and techniques for densities ranging from sub-nuclear to 
, a few times saturation density. Accordingly, and spurred by the study of rare isotopes, future theoretical research is 
' also directed towards a systematic study of neutron-rich matter. 

A unique method to probe the internal composition of a neutron star is by tracing its temperature evolution with 
' cooling simulations, see e.g., [3-6]. Neutron stars are born in supernova explosions and the interior temperatures 
^f-^ initially exceed T ~ 10^^ K. As the neutron star cools, neutrinos begin to free stream and essentially leave the star 
• without further interaction. As a consequence, after about 30 s, the long-term cooling of neutron stars is controlled by 
\ neutrino emission. This stage lasts up to about 10^ years of age, when cooling by emission of photons becomes more 
. effective. Importantly, the neutron stars remain luminous enough during the cooling, so that the surface temperature 
,_H ' can be extracted from space telescope data and the theoretical cooling curves can be confronted with observations. 
^ , For proton fractions np/(n„ -|- Up) < 1/9, direct beta decay does not proceed in neutron star matter due to the 
imbalance of the neutron and proton Fermi momenta. Therefore, and in the absence of nucleon superfluidity, the 
dominant neutrino emission comes from vi? bremsstrahlung in nucleon-nucleon collisions and the so-called modified 
Urea process. The latter corresponds to in-medium beta decay, where the momentum difference between the decaying 
neutron and the final proton is absorbed by scattering off a second nucleon. For low temperatures T ^ Ep^ (where 
Ep denotes the Fermi energy and we use fc^ = 1), the temperature dependence of these different emission channels 
is easily power-counted via the degeneracy of the fermions involved in the emission. It follows that the emissivity e, 
which is the total neutrino energy emitted per unit volume and unit time, scales as e ^ . 

However, the density dependence and strength of the emissivity requires, as input, accurate calculations of the 
nucleon-nucleon scattering amplitude in the many-body medium. Specifically, the bremsstrahlung processes involve 
the scattering amplitude for nucleons on the Fermi surface, whereas the modified Urea process involves also off- 
shell scattering, when e.g., the scattered neutron beta decays and thus its momentum is approximately given by the 
proton Fermi momentum. A further theoretical challenge lies in the fact that the dominant contribution to neutrino 
emissivities comes (for bremsstrahlung in nn collisions exclusively) from the non-central parts of the nucleon-nucleon 
amplitude, in particular from the tensor force [7]. In addition to neutrino emissivities, the general importance of 
non-central interactions for neutron star properties has been revived recently, where the effects on the magnetic 
susceptibility [8] as well as on P-wave pairing in neutron star cores [9] have been demonstrated to be crucial. 

For the nucleon scattering amplitude, the benchmark calculation of Friman and Maxwell takes into account the 
long-range one-pion exchange tensor force explicitly and estimates the effects of short-range correlations by cutting off 
the interaction at short distances and including the short-range rho-exchange tensor force [7]. In this work, we start 
from the free-space low-momentum nucleon-nucleon interaction Viowfe [10,11]. The construction of Viowfe is motivated 
by the differences of the realistic nucleon-nucleon potential models at short distances. The model dependence at short 
distances r < d « 0.5 fm originates from the fact that the interaction cannot be resolved from scattering experiments 
probing low momenta p < A w 2.0 fm^^ (where A = l/c?)- Note that the realistic potential models are fitted to phase 
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shifts below i?iab ~ 350 McV, corresponding to A « 2.1 fm~ . A systematic method to remove the model dependence 
is provided by the renormalization group (RG), where the high momentum modes with p A are integrated out to 
construct the physically equivalent effective theory. The renormalization group in this context is used as a tool to 
guarantee that the phase shifts are preserved by the low-momentum interaction under the renormalization. Since 
Viowk does not have momentum components larger than the cutoff A, it does not have a strongly repulsive core. As 
a consequence, one does not have to compute a Brueckncr G matrix from Viow k in many-body applications. In the 
sense of the RG, short-range correlation effects are implicitly included in Viowfe [12]. 

In addition to the phenomenological short-range correlation effects discussed by Friman and Maxwell, the scattering 
amplitude is screened by the particle-hole polarization of the many-body medium. For the in-medium tensor force, 
particle- hole screening effects are very important. This follows from a general spin-recoupling argument due to the 
interference of the central spin-spin part and the tensor part of the nucleon-nucleon interaction, leading to a substantial 
decrease of the latter [9,12]. Moreover, the presence of the Fermi sea defines a preferred frame, which leads to novel 
non-central parts in the effective interaction and the scattering amplitude in the many-body medium [9,13]. In pure 
neutron matter, the scattering amplitude on the Fermi surface has been computed to second-order in Viowfe, with 
particular attention to the spin-dependence and non-central interactions, where it was found that the particle-hole 
screening leads to a substantial decrease of the tensor force and a significant long-wavelength center-of-mass tensor 
force induced by the medium [9] . As our understanding of the renormalization of the nucleon-nucleon interaction in 
dense matter improves, it is thus important to include the in-medium modifications of the nuclear force as well as 
the novel non-central contributions in the calculation of the neutrino emissivities. Furthermore, with recent results 
for the neutron superfluid S- and P-wave pairing gaps including polarization effects obtained in [9,14], a consistent 
calculation of the emissivity is needed. 

The neutrino emissivities also receive contributions from the spin-orbit force, which were included in recent cal- 
culations of Hanhart et al. [15] and van Dalen et al. [16] starting from the free-space scattering amplitude. As for 
the tensor force, particle-hole polarization effects reduce the spin-orbit interaction in neutron matter [9]. As a con- 
sequence, the superfluid pairing gaps in neutron star cores are strongly suppressed to below few keV at second 
order in Mow ft for the pairing interaction [9]. Therefore, one expects that the transition to the P-wave superfluid 
phase of neutrons is only reached in the very late stages of neutron star cooling. The consistency with data in present 
cooling simulations of Yakovlev et al. also requires low critical temperatures of the '^P2 superfluid, Tc < 2 • 10^ K [6]. 
This corresponds to an angle-averaged gap (as calculated in [9]) in the ni j = state of A < 30keV. We note that on 
the level of the free-space nucleon-nucleon interaction considerably larger P-wave pairing gaps of A w 0.35 MeV are 
predicted at nuclear matter density [17]. 

If we thus consider neutron star matter above nuclear matter density, the ^Sq superfluidity of neutrons ceases to 
exist due to the repulsion in the nuclear force and one expects that, at the relevant core temperatures, the neutrons 
are in the normal phase and the protons are superconducting. In this case, the modified Urea process is strongly 
suppressed due to proton superfluidity, and the dominant neutrino emission process will be neutrino bremsstrahlung 
in nn collisions. This is the motivation to focus on the in-medium modiflcation of the bremsstrahlung rate in this 
work. The second motivation comes from the fact that even at lower densities, where both neutrons and protons pair 
in the ^ So channel, the neutrino pair emissivity is more effective compared to the modified Urea process [18] , although 
this conclusion depends on the values of the pairing gaps and neutrino emission mainly proceeds through the so-called 
Cooper Pair-Breaking and Formation (PBF) process in the regime 0.2 Tc < T < Tg [19]. 

In this work, we use as input the second-order results for the neutron-neutron scattering amplitude on the Fermi 
surface computed in [9]. We start by giving the general spin-dependence of the scattering amplitude on the Fermi 
surface, with a short discussion of the novel contributions in Section II. The non-central interactions are included 
in the calculation of the emissivity from neutrino bremsstrahlung in Section III. This Section follows closely the 
derivation of Friman and Maxwell [7] . The results for the emissivity over a range of densities is given at the end of 
Section III. Finally, we conclude with a discussion of superfluidity in neutron stars, which takes into account recent 
results for the pairing gaps [9,14]. We present some general arguments for the constraints on the density dependence 
of the gaps and compare the bremsstrahlung rate to the PBF process for various temperatures. Revised estimates of 
the np and pp bremsstrahlung and modifled Urea rates will be reported in a subsequent publication, as they require 
the nucleon-nucleon scattering amplitude off the Fermi surface as well as an extension to asymmetric matter. 

II. SPIN-DEPENDENCE OF THE SCATTERING AMPLITUDE IN NEUTRON MATTER 

The scattering amplitude for neutrons on the Fermi surface contains the free-space central (scalar and spin-spin) 
and non-central (spin-orbit and tensor) parts. In addition, the many-body medium can induce effective interactions, 
which depend on the two-body center of mass momentum P = pi -\- p2 = pz + Pi and one has in general [9] 



2 



A„,,„M, q', ^) = E ?') OU^„,{q, q', P) (1) 

i 

= Acalar(<?, (?') + ApinCg, CTi • CTa + Apin-orbit (?') «(cri + CT2) ■ q y. cj' 

+ Atcnsor{q,q') Si2{q) +AeKch. tensor (g, (?') '5'i2(g') + Acm tensor (g, S'i2(P) 
vector 

(q, q') i{(Ti - (T2) ■ q X P + Ac ross vector 

{q,q') {(Ti X CT2) ■ (q' X P), (2) 

where q = Pi — Ps and q' = Pi — P4 denote the momentum transfers in the direct and exchange channel respectively, 
the tensor operator is defined as Si2{q) = <ti • q <T2 • q — 1/3 cti • cr2, and the scattering amplitude is defined in units of 

the density of states to* fci?^/7r^. For particles on the Fermi surface, the momentum transfers and the ccnter-of-mass 
momentum are orthogonal and one has + q'"^ + = 4fc^^. Therefore, the various parts Ai depend only on the 
magnitude of the momentum transfers ^=191 and q' = \q'\. We also note that the non-central operators are defined 
with unit vectors [8]. Finally, the tensor operators given in Eq. (2) are linearly dependent, with 6*12(9) + Suiq') + 
Si2{P) = 0. In the second-order calculation all tensors are kept explicitly [9], and for the emissivities we then eliminate 

5'l2(P), leading to Acnsor = Acnsor - Am tensor aS WCU aS Axch. tensor = Axch. tensor - Am tensor- In the Scattering 

amplitude, both direct and exchange terms are accounted for; e.g., to lowest order, one has A = V\owk — Po-PfeMow/t, 
with spin- and momentum-exchange operators Per and Pfe. The latter two operators in Eq. (2) do not conserve the 
spin of the interacting particle pair and are induced in the medium due to the screening by particle-hole excitations. 

Although it is generally known that particle-hole polarization effects arc very important in nuclear physics, this is 
the first calculation where these are included for neutrino emissivities. In the following we compute the emissivity 
without superfluid effects, in order to assess the renormalization of the tensor and spin-orbit forces in matter and 
the contributions of the novel non-central forces, where the cm tensor is included in the conventional tensor parts. 
Subsequently, we discuss polarization effects for the superfluid properties and compare the bremsstrahlung rate to the 
dominant PBF process [19]. In both parts we give results for the density dependence of the emissivity. 




FIG. 1. Feynman diagrams contributing to the neutrino emissivity from bremsstrahlung to second-order in Viowfc- The Viowfc 
vertex includes both the direct and the exchange term and the dashed line corresponds to the neutral current. As in [7], the 
emissivity includes permutations of the neutral current attached to all external lines, whereas the coupling to the internal 
nucleon lines is suppressed, because it does not lead to the small energy denominator in the additional nucleon propagator as 
in Eq. (5). 



III. NEUTRINO BREMSSTRAHLUNG IN NEUTRON-NEUTRON COLLISIONS 

Following Friman and Maxwell [7] , the emissivity from neutrino pair bremsstrahlung in neutron- neutron scattering 
is given by (for h = c = 1) 

where Pi denote the momenta of the incoming and outgoing neutrons and Qi^2 = ('^1,2, Qi,2) label the neutrino energies 
and momenta. The delta functions account for energy and momentum conservation, and = + UI2 is the total 
neutrino energy. N^, denotes the number of neutrino species and s = 2 is a symmetry factor for the initial neutrons, 
when the emission occurs in the final state or vice versa. The function J^{Ep.) = f{Ep^) f(Ep^) (1 — f{Ep^)) (1 — 
f{Ep^)) is the product of Fermi-Dirac distribution functions f{E) = {ex.p{E/T) + 1)~^, with neutron energies Ep^. 
The matrix element A^„„ includes the nucleon- nucleon scattering part and the coupling to the emitted neutrino pair. 
For the bremsstrahlung process the corresponding Feynman diagrams at tree-level and with second-order contributions 
are shown in Fig. 1. The second-order particle- hole intermediate states include all possible excitations for interacting 
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particles on the Fermi surface. In the particle-particle channel, the cutoff in Viowfc provides a regulator, and we 
evaluate the phase space (including hole-hole states) exactly without angle-averaging approximation. We note that 
the logarithmic contribution to the quasiparticle interaction (due to the BCS singularity) is integrable and leads to 
finite Fermi liquid parameters [20], and thus finite emissivities. 

Neutrino pair emission from a neutron line is given by the neutral current V — A weak interaction 

-^neutral = " ^ Xl (^/^O - 9 A S^,i ai) X2 I,, , (4) 

with Fermi coupling constant Gp — f .166 • 10^'"' GeV~^, the neutrino current = u{Qi) 7^ (1 — 75) u{Q2) and the 
weak axial-vector coupling constant gA = 1.26. The non-relativistic nucleon spinors are denoted by x and m(Qi,2) are 
relativistic spinors for the neutrinos, which are taken to be massless. As in [7] , we use a non-relativistic approximation 
for all nucleon propagators, where the lowest term in an expansion in inverse powers of the nucleon mass is retained. 
In addition, one neglects the neutrino pair energy uj^ compared to the Fermi energy, since the emitted neutrinos are 
thermal. For the nucleon propagator G, this approximation yields 

iG{p±Q,,Ep±w^)=±iu)-\ (5) 

where the positive sign holds if the weak current is attached to an outgoing nucleon, negative otherwise. It follows 
that, in the non-relativistic approximation for the nucleon propagators, the vectorial part of the neutral current does 
not contribute [7]. 

The matrix element A4nn includes the strong interaction part with spin-dependence given by Eq. (2). The spin 
sum over the squared matrix element is carried out independently for the pieces coming from the nucleon-nucleon 
amplitude and the weak interaction, where only the non-central parts in the amplitude are found to contribute. After 
contraction with the lepton trace given by 

Tt(/i it) = 8 (Qi, Q2, + Q2, Qi, - Qij Qi-Q2+i Ciaj^ QI) , (6) 

one finds 

^4, 



spin 



= 64 5^ G\ 2 IPiUT' ( Aensor(«) S') + -^exch. tensor(9. «') " Aensor(5, «') Axch. tensor(5, (l) 

+ -^spin-orbit (9; 1 ) + •^diff. vector 9 ) + 3 .Aj-j-Qsg vector (9) 9)1) (7) 



where we have dropped terms that vanish upon angular integrations over Qi,2 in the emissivity, Eq. (3), when the 
neutrino momenta are neglected in the momentum-conserving delta function. 

For the evaluation of the phase space, one can easily perform the integrals over the neutrino momenta with |Qi,2| = 
wi,2 by inserting 



1 = y 8{lOv - wi - W2), (8) 

as the energy-conserving delta function depends only on the total neutrino energy. Moreover, one can decouple the 
angular parts in the neutron phase space and trade the radial momentum for energy integrals, by restricting the 
interacting neutrons to the Fermi surface, since they are strongly degenerate for typical neutron star temperatures. 
Corrections to this approximation scale as Tj Ef„ ■ For this purpose, one replaces 



d-'pi ^ rfVi Sip, - kp^ ) J dEp^ . (9) 

Finally, we introduce the integration over momentum transfers through respective delta functions in the direct and 

the exchange channels, 

1 = j d''qS^q-pi+p3) = j dV'5^(g'-Pi+P4). (10) 
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After carrying out the angular integrations, we find the general expression for the emissivity 



where I^p is the result of the integrals over the total neutrino and neutron energies convoluted with the thermal 
distribution factors J-'{Ep.) which is identical to the expression in [7], and the square of the scattering amplitude 
averaged over the Fermi surface is given by 

dq f dq' fcF„e(4fc|, - q^ - q'^) ,2 , f dq 



"0 " \ F,, <i V ) " 

The remaining two-dimensional integral is calculated numerically and our results for Ni, — 3 neutrino flavors are shown 
in Fig. 2 for densities ranging fromfci. =1.0-2.0fm~\ All results given in Fig. 2 include the effective mass obtained 
from the lowest order Vlowfe- The effective mass varies from to*/to„ = 0.95 at kp^ — 1.0 fm^^ to Trt^/mn = 0.78 at 
kpn = 2.0 fm~^, and is in this range well approximated by a linear curve versus Fermi momentum. We note that 
one expects an increase of the effective mass in the induced interaction (which is compensated by the quasiparticle 
strength 2:^.^), as can be seen from the results of the full RG calculation for neutron matter [14]. Furthermore, our 
results do not include a renormalization of qa in the medium to w 1.0. 
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FIG. 2. The neutrino emissivity from bremsstrahlung in neutron-neutron scattering e„„ versus Fermi momentum kp^ in 
neutron matter. All curves include the lowest-order effective mass obtained from l^owfc, see [9]. The curve labeled V denotes 
the lowest-order emissivity obtained from the free-space low-momentum interaction. In comparison, we give the results obtained 
by Friman and Maxwell (curve F) using the uncorrelated one-pion exchange (OPE) tensor force without the exchange terms [7]. 
(In [7], the inclusion of exchange terms, the p tensor force, and short-range correlation effects yields a multiplicative suppression 
factor of « 0.62 at nuclear matter density or kp^ = 1.7fm~^.) The thin curves correspond to the Viowfc result with, respectively, 
second-order renormalization of the tensor, spin-orbit or spin non-conserving forces in the medium included. The curve labeled 
T is the full second-order result, where both particle-hole channels and the particle-particle channel are taken into account. 



First, we compare the lowest-order Viowfc results to the calculation of Friman and Maxwell from uncorrelated one- 
pion exchange (OPE), see Eq. (52) in [7]. We remark that exchange terms, the inclusion of the p tensor force and 
correlation effects lead to a multiplicative suppression factor in the calculation of Friman and Maxwell of k, 0.62 
at nuclear matter density. We find that T^ow fc gives similar rates, but without the need to estimate the correlation 
distance, which is experimentally unconstrained in neutron matter. 

Recently, Hanhart et al. [15] and van Dalen et al. [16] also computed the neutrino pair emissivity from 
bremsstrahlung employing Low's theorem for soft emission, with the free-space on-shell T matrix as input. These 
results provide a model-independent low-density limit on the emissivity from neutrino bremsstrahlung. In both works, 
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the cmissivity was found to be reduced by a factor w 1/4 compared to the OPE result of Friman and Maxwell at 
saturation density. However, the applicability to relevant neutron star densities is limited when one starts from the 
free-space scattering amplitude. This is due to the fact that at second-order in the T matrix the spin-spin and tensor 
part of the amplitude mix due to screening in the particle-hole channel. If one denotes the tensor part of the free- 
space amplitude by Ttonsor, then the second-order contribution will be proportional to Ttensor fc_F„ «s at low momenta, 
where as is the S-wave scattering length c;oming from the spin-spin part of the T matrix. Even at low-densities, e.g., 
fcf^ = 1/2 fm~^ (i.e., p = 1/40 po), this would be a very large correction, which is not accounted for in [15,16]. 

Beyond the lowest-order result, we find that the renormalization of the non-central parts of the nucleon-nucleon 
interaction in the medium considerably reduces the cmissivity, especially at sub-nuclear densities, as shown in Fig. 2. 
We emphasize that a second-order calculation cannot give final results, but it provides a range for the effects due to 
particle-hole screening. Prom the denominators of the induced interaction [22] , one expects higher orders to somewhat 
decrease the second-order results. However, this argument is not straightforward for non-central interactions and a 
definitive conclusion requires explicit calculations of higher-orders with full non-central spin-dependence. Referring 
to a previous calculation of the induced quasiparticle interaction (i.e., for q = Q) including tensor forces [23], it was 
also found that the tensor force is reduced in the medium. For the cooling of neutron stars, as well as spin-isospin 
response in supernovae, the modification of transport properties in the medium is very important. Our work shows 
that particle-hole effects must be taken into account in a realistic calculation of the emissivities. 

As discussed in the Introduction, an accurate estimate of the bremsstrahlung rate is important particularly at high 
densities, since the competing processes except for possible emission by PBF are expected to be strongly suppressed 
due to proton superfluidity. Therefore, we will compare the strength of the open PBF channel to the emission from 
a non-superfluid core of neutrons in the next Section. 

IV. SUPERFLUIDITY AND COMPARISON WITH PBF PROCESSES 

The temperatures in the interior of cooling neutron stars can be well below the critical temperature for neutron or 
proton superfluidity. This leads to a strong reduction of the neutrino emissivities since the fraction of particles that 
are unpaired scales exponentially with the temperature as exp(— 2A/r), where A denotes the zero temperature gap. 
Thus, the modification of the emissivities relies on an accurate calculation of the proton and neutron pairing gaps in 
neutron star matter, which must include the particle-hole polarization in the medium. 

Before discussing recent results for the neutron pairing gaps in pure neutron matter [9,14], we proceed with some 
remarks on the superfluid properties and rather general in-medium modification arguments. The strongest attraction 
in the nuclear force is in the S-wave for laboratory energies below E'lab ^ 250 — 260 MeV, i.e., for back-to-back 
scattering of particles with momenta k = kp^ ^ 1.7 — 1.8 fm"'^. Thus, for densities below nuclear matter density 
P ^ Po, one expects that both neutrons and protons form a superfiuid in the isotriplet ^Sq channel (Note that for 
typical proton fractions of np/(n„ -|- rip) « 0.05, one has for the proton Fermi momentum kp^ « 1/3 kp^)- At higher 
densities, one concludes from the free-space scattering phase shifts that the neutrons are expected to pair in the ■^P2 
state [21]. Eventually, the ground state of matter at high densities has to be determined in a model of dense matter, 
as realistic interactions are constrained to relative momenta A; < 2.1 fm~^. 

It is well known that polarization effects on the nucleon-nucleon interaction, which arc necessary in order to satisfy 
the Pauli principle in microscopic calculations [22], lead to a strong reduction of the superfluid gaps in neutron 
stars [9,14,24-26]. The effect of the induced interaction on pairing can be understood by considering the second-order 
contributions. Higher-order terms modify the strength of the second-order result, but usually do not alter whether 
the induced interaction is attractive or repulsive in the particular channel. 

For S-wave pairing of neutrons, the dominant part in the quasiparticle interaction is the central spin-spin delta 
function Go « 0.6 — 0.8 [14]. Neglecting the smaller contributions, one has for the induced pairing interaction in the 
5 = state, 

-4j?ontral(cOS(?,) = 3 {Uiq/kpJ + Uiq'/kpj), (13) 

where direct and exchange particle-hole channels are accounted for, U{q/kp^) denotes the (positive) static Lindhard 
function and for back- to-back scattering q = kp^ yj2 — 2 cos^q {q' = kp^^2 + 2cos6q) with scattering angle 9q. The 
projection of the sum of Lindhard functions in Eq. (13) on S-wave yields {U{q/kp^) + U{q' /kp^))i=o = 2(l + 21og2)/3 w 
1.59 (see also [27]). As a result, spin fluctuations will reduce the pairing interaction and consequently the superfluid 
gap will close at somewhat lower densities than one would expect from the free-space ^Sq phase shifts. In the RG 
calculation of the effective interaction and the scattering amplitude on the Fermi surface [14] (for a discussion of the 
RG approach see also [28]), it is found that the maximum ^Sq pairing gap is reduced to 0.8 MeV at kp^ ~ 0.8 fm~^ 
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and that the gap disappears at fc^?^ « 1.5 fm or p « 2/3po- This corresponds to a maximum critical temperature 
of = 5.3- lO^K. 

The study of polarization effects on P-wave pairing at higher densities is more involved, since spin-orbit and tensor 
forces are crucial for the reproduction of the P-wave phase shifts in vacuum and consequently one has to address 
the renormalization of non-central interactions in the medium. In particular, the medium-induced spin-orbit force 
leads to a strong suppression of the ^P2 gaps, which is e.g., due to the interference of the central spin-spin and the 
spin-orbit force at second order [9] . A similar argument as for S-wave pairing reproduces the effect of induced pairing 
interactions qualitatively and one has 

4plorbit(cos0g) = Go {Vso) {U{q/kFj + U{q'/kF„)) g^, (14) 

where the dimensionless (Vso) < denotes an averaged spin-orbit matrix element. The induced contributions due 
to the mixing of spin-orbit and tensor forces are also repulsive, with a similar but more complicated momentum 
dependence. For the induced '^P2 pairing matrix element, the largest contribution comes from the same I ~ 
projection of the sum of the Lindhard functions, where the additional momentum-dependent factors qq' are absorbed 
in the spin-orbit operator L ■ S ^ i{cri + ■ q x q' (note the unit vectors in Eq. (2)). In [9] it was found that 
second-order polarization effects lead to a strong depletion of the ^P2 pairing gap from A w 0.3 MeV (obtained from 
the free-space Viowfe) to superfluid gaps on the level of few keV at nuclear matter density. This corresponds to a 
ratio of 0.45 of the second-order to lowest-order pairing interaction, and therefore it is expected that the suppression 
of the gap at second order is robust with pairing gaps below A < 1 — lOkeV. The latter value corresponds to a 
critical temperature of Tc — 10^-*"^-*K. With surface temperatures from observational data ^ lO'^-^K, i.e., core 
temperatures T ^ 10'''^K (see e.g.. Table 2 in [29]), it is thus expected that the ^P2 superfluid phase is only reached 
at late cooling stages. In fact, Yakovlev et al. have checked that for critical temperatures Tc < 2 • 10* K the ^P2 phase 
has no impact for the cooling of middle-aged neutron stars [6] . We therefore proceed and compare the bremsstrahlung 
emissivity in normal matter to the emission of neutrinos from the ^Sq superfluid condensate only. 
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FIG. 3. Comparison of the emissivity from neutrino pair bremsstrahlung in non-superfluid neutron matter (total in Fig. 2) 
and the ^So PBF process as a function of Fermi momentum and for various temperatures. Results are shown using the ^So 
superfluid gaps obtained in the RG approach with adaptive Zkp factor [9]. 



Although superfluidity strongly suppresses the standard neutrino emission channels, there is a powerful mechanism 
for neutrino emission in superfluid matter due to the PBF process [19]. This is a significant source for neutrinos 
from the density range, where the temperature lies between Q.2Tc{kF,^) < Tc{kp^). The emissivity from the ^Sq 
neutron pair-breaking and formation process is given by 

ePBF = 1.170 • 10^1 TV, ^ ^ Tj F(t) erg cm'^ s'^ (15) 
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where the function F{t) depends on the critical temperature r = T/Tc{kp^) at given neutron Fermi momentum (for 
S-wave pairing Tc{kp^) ~ 0.57A(fci7'^)). Details and a parametrization of F{t) which wc employ can be found in [18]. 
In Fig. 3 we compare our results for the bremsstrahlung emissivity in the normal state to the ^So PBF process. While 
it demonstrates that the PBF process is extremely effective at higher temperatures T > 10*'^ K, at lower temperatures 
the density range for emission via PBF is rather narrow and we expect the volume-integrated bremsstrahlung emission 
to dominate. Fig. 3 also nicely illustrates how the cooling of neutron stars is able to probe the internal structure, 
with the enhanced emissivity from the PBF process serving as a clear signal of the superfluid phase. 

We finally note that at higher neutron densities the proton PBF process is also in effect. However, at lower 
temperatures the cooling through neutrino bremsstrahlung in normal-state nn collisions will be larger for the same 
reason as above, and in addition due to the smaller number of protons. 

V. SUMMARY 

Wc have computed the emissivity from neutrino pair bremsstrahlung in nuclcon-nucleon collisions in pure neutron 
matter, within an effective theory of quasiparticle interactions in the vicinity of the Fermi surface. The effective 
scattering amplitude is calculated from the model-independent low- momentum nucleon-nucleon interaction VJowfc to 
second order, keeping the full non-central spin dependence [9]. We find that inclusion of medium modifications, in 
particular the renormalization of the tensor force, reduces the emissivity compared to the tree-level Viowfc result by a 
multiplicative factor 0.64 at nuclear matter density (or a factor 0.5 relative to the direct one-pion exchange estimate). 
At sub-nuclear densities, the reduction is 0.2 at a fifth of nuclear matter density. Furthermore, we find that the effect 
of spin non-conserving parts in the scattering amplitude is rather small. While the temperature dependence of the 
emissivity is naturally very important for cooling simulations, the density dependence needs to be considered as well, 
since the luminosity of the neutron star involves an integral over the volume of the star. 

When polarization effects are included in the ^P2 neutron pairing interaction, a considerable reduction of the 
superfluid gaps was found [9]. Therefore, the neutrino bremsstrahhmg process could be more important for the cooling 
of neutron stars than believed, since superfluidity of protons suppresses Urea processes as well as bremsstrahlung in np 
and pp scattering, whereas neutrons remain in the non-superfluid phase at higher densities at typical core temperatures 
in the neutrino cooling stage. Even at lower densities, where neutrons form a ^So superfluid, the lore is that neutrino 
bremsstrahlung from neutrons dominates the modified Urea process in the presence of superfluidity [18], although 
the PBF channel is considerably more effective. We also note that the bremsstrahlung process, in contrast to the 
Urea channels, produces u^j, and Vt neutrinos. In order to compare our results for the bremsstrahlung rates to the 
emission through the PBF process, we have shown results for the PBF process for realistic pairing gaps obtained in 
the RG approach [14]. This nicely demonstrates that the PBF channel is very effective and can act as a powerful 
signal of superfluidity for temperatures comparable to the maximal critical temperature, but is restricted only to a 
narrow density range for lower temperatures. In the latter regime, one thus expects the integrated emissivity from 
bremsstrahhmg to dominate. 

Our results provide neutrino bremsstrahlung rates derived from successfully used effective nuclear interactions. The 
emissivities can be used as microscopic input for neutron star cooling simulations, in conjunction with the superfluid 
gaps calculated in [9,14]. A RG calculation of effective interactions in asymmetric matter, which will address both 
the effects of induced interactions on proton pairing in neutron star matter as well as higher-order contributions to 
non-central interactions, is in preparation [30]. A self-consistent treatment of the tensor force is necessitated by a 
substantial renormalization at second order. Such studies will further constrain neutrino emissivities microscopically 
in a consistent framework. Our results can then be incorporated in cooling simulations of neutron stars to possibly 
constrain the structure of neutron stars and their densest interiors. 



ACKNOWLEDGMENTS 

We thank Bengt Priman and Dick Furnstahl for useful discussions. AS is supported by the NSF under Grant No. 

PHY-0098645 and through an Ohio State University Postdoctoral Fellowship. The work of PJ and CG is supported 
jointly by the Natural Sciences and Engineering Research Council of Canada and the Fonds Nature et Technologies 
of Quebec. 



8 



[1] C.J. Pethick and D.G. Ravcnhall, Ann. Rev. Nucl. Part. Sci. 45 (1995) 429. 

[2] H. Heiselborg and V.R. Pandharipandc, Ann. Rev. Nucl. Part. Sci. 50 (2000) 481. 

[3] C.J. Pethick, Rev. Mod. Phys. 64 (1992) 1133. 

[4] D. Page, in "The Many Faces of Neutron Stars", Proceedings of the NATO Advanced Study Institute, Eds. R. Buccheri, 

J. van Paradijs, M.A. Alpar, Kluwer, 1998, astro-ph/9706259. 
[5] A.D. Karninkcr, D.G. Yakovlcv and O.Y. Gnedin, Astron. Astrophys. 383 (2002) 1076. 
[6] D.G. Yakovlev, O.Y. Gnedin, A.D. Kaminker, K.P. Levenfish and A.Y. Potekhin, astro-ph/0306143. 
[7] B. Friman and O.V. Maxwell, Astrophys. J. 232 (1979) 541. 
[8] E. Olsson and C.J. Pethick, Phys. Rev. C66 (2002) 065803. 

[9] A. Schwenk and B. Friman, accepted for publication in Phys. Rev. Lett., nucl-th/0307089. 
[10] S.K. Bogncr, T.T.S. Kuo, A. Schwenk, D.R. Entcm and R. Machleidt, Phys. Lett. B576 (2003) 265. 
[11] S.K. Bogner, T.T.S. Kuo and A. Schwenk, Phys. Rep. 386 (2003) 1. 
[12] A. Schwenk, G.E. Brown and B. Friman, Nucl. Phys. A703 (2002) 745. 
[13] J.L. Forest, V.R. Pandharipande and J.L. Friar, Phys. Rev. C52 (1995) 568. 
[14] A. Schwenk, B. Friman and G.E. Brown, Nucl. Phys. A713 (2003) 191. 
[15] C. Hanhart, D.R. Philips and S. Roddy, Phys. Lett. B499 (2001) 9. 
[16] E.N.E. van Dalen, A.E.L. Dieperink and J.A. Tjon, Phys. Rev. C67 (2003) 065807. 
[17] M. Baldo, 0. Elgar0y, L. Engvik, M. Hjorth-Jensen and H.-J. Schulze, Phys. Rev. C58 (1998) 1921. 
[18] D.G. Yakovlev, A.D. Kaminker, O.Y. Gnedin and P. Hacnscl, Phys. Rep. 354 (2001) 1. 
[19] E. Flowers, M. Ruderman and P. Sutherland, Astrophys. J. 205 (1976) 541. 
[20] O. Sjoberg, Nucl. Phys. A209 (1973) 363. 

[21] M. Hoffberg, A.E. Glassgold, R.W. Richardson and M. Ruderman, Phys. Rev. Lett. 24 (1970) 775. 
[22] S. Babu and G.E. Brown, Ann. Phys. 78 (1973) 1. 

[23] W.H. Dickhoff, A. Facsslor, H. Muothcr and Shi-Shu Wu, Nucl. Phys. A405 (1983) 534. 

[24] J.W. Clark, C.-G. KaUman, C.-H. Yang and D.A. Chakkalakal, Phys. Lett. B61 (1976) 331. 

[25] J. Wambach, T.L. Ainsworth and D. Pines, Nucl. Phys. A555 (1993) 128. 

[26] C.J. Pethick and D.G. Ravenhall, Ann. N.Y. Acad. Sci. 647 (1991) 503. 

[27] H. Heiselberg, C.J. Pethick, H. Smith and L. Viverit, Phys. Rev. Lett. 85 (2000) 2418. 

[28] A. Schwenk, B. Friman, G.E. Brown, in "Hirschcgg 2003, Nuclear Structure and Dynamics at the Limits", nucl-th/0302081. 
[29] D.G. Yakovlev, K.P. Levenfish and Yu.A. Shibanov, Physics Uspekhi 42 (1999) 737. 
[30] A. Schwenk, B. Friman and R.J. Furnstahl, in preparation. 



9 



